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Abstract 

We study the scattering of waves off a potential step in deformed honeycomb lattices. For small 
deformations below a critical value, perfect Klein tunneling is obtained. This means that a potential 
step in any direction transmits waves at normal incidence with unit transmission probability, 
irrespective of the details of the potential. Beyond the critical deformation, a gap in the spectrum 
is formed, and a potential step in the deformation direction reflects all normal-incidence waves, 
exhibiting a dramatic transition form unit transmission to total reflection. These phenomena are 
generic to honeycomb lattice systems, and apply to electromagnetic waves in photonic lattices, 
quasi-particles in graphene, cold atoms in optical lattices. 
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Scattering of relativistic fermions is fundamentally different from that of non-relativistic 
ones, since relativistic fermions are described by the Dirac equation which is first order 
in momentum (rather than the second order Schrodinger equation). Relativistic fermions 
(massive or massless) incident normally upon a potential step of height Vq, exhibit non-zero 
transmission probability, even when their energy is smaller than Vq [1] . This behavior of rel- 
ativistic fermions is called Klein tunneling and stands in sharp contrast to the more intuitive 
result of quantum mechanics for non-relativistic particles, where the transmission probabil- 
ity vanishes completely when the height of the step is greater than the particle's energy. 
Outside the step, the state is of positive energy, whereas inside the step there is a propa- 
gating negative energy state. This unique scattering process, has never been experimentally 
verified, since an experiment designed to observe Klein tunneling with elementary particles 
requires high fields which are not currently available. However, it has been predicted that 
the charge carriers in graphene, that obey the massless Dirac's equation, exhibit similar 
behavior [2]. More specifically, it has been suggested that charge carriers in graphene expe- 
rience non-resonant unit transmission in monolayer graphene, and total reflection in bilayer 
graphene [3]. Experiments with bipolar junctions were able to show very high conductance 
in the presence of a gate voltage, indicating high transmission probability [4-6] . In addition, 
unusual transmission properties are predicted in honeycomb photonic crystals [7, 8]. 

Here we study the dynamics of waves in a deformed honeycomb photonic lattice, and in 
particular the tunneling process into a refractive index step. We find that, up to the critical 
deformation in which a gap in the spectrum is formed, non-resonant unit transmission is 
obtained at normal incidence, irrespective of the details of the potential. That is, surpris- 
ingly, deformed honeycomb lattices also display perfect tunneling as well as non-deformed 
honeycomb lattices. At deformations stronger than the critical one, we find non-resonant 
total reflection in the deformation direction. Generally, optical structures exhibiting unit- 
transmission are resonant, and are characterized by transmission peaks, e.g., Fabry-Perot 
etalon. Thus, our study introduces a new domain of light transport in photonic struc- 
tures, displaying non-resonant effects, and offering an opportunity to directly observe Klein 
tunneling. 

Paraxial propagation of a monochromatic field envelope with a honeycomb refractive 
index is given by [9] 

■ di/j 1 ^2 , k6n{x,y) , , 

^aJ = -2^^-^-^^^ = ^^' 
where 6n is the modulation in the refractive index, k is the wave-number and uq the back- 
ground refractive index. Note that k6n{x, y)/no is an optical potential with an opposite sign, 
which means that light is attracted to higher index of refraction. Since 6n is independent 
of z, Eq.(l) has solutions of the form ip^x, y, z) = U{x, y) exp(z/?z), where (3 is the propaga- 
tion constant, and f/ is a solution of HU = (3U, hence (3 is analogous to the energy. Like 
other periodic systems, the system can be described by tight binding (TB) and hence the 
analogy with graphene. The honeycomb lattice is comprised of two triangular sub-lattices. 
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denoted by A and B (Fig. la). Therefore, when writing the Hamiltonian in Wannier basis, 
one must consider two sets of amphtudes, an, bn, associated with the two sites A and B in 
each unit cell located at n. Assuming only nearest neighbors hopping (coupling), the TB 
Hamiltonian resulting from Eq.(l) reads [10] 

Ho = -Yl {<bri+s, + bn+s^a^n) , (2) 

where t/s are the hopping parameters, and 5j are the vectors connecting the nearest neigh- 
bors. The anisotropy is manifested by different hopping parameters in different directions, 
and can be realized by various means [11-13]. We consider a photonic lattice with different 
spacing between the sites (Fig. lb). Such deformations are uniaxial, i.e., t2 = = t and 
ti = '-ft, where 7 > 1. We refer to the x— direction as the deformation direction. Expending 
ttn, bn in Fourier space and defining ipk = J2j exp{idjk), we can write 

Ho = ^Yl = + sXm{v^fc}, (3) 

BZ 

where cTj are Pauli matrices and = (a]^ b^,) is a pseudospinor. The spectrum is obtained 
by the eigenvalues of Hk yielding positive and negative branches: 

f3 = ±t\l 2 + 72 + 47 cos ^ cos + 2 cos"fc^, (4) 

where a is the lattice constant. For 7 < 2, the two branches intersect at two inequivalent 
points in the first Brillouin zone, known as the Dirac points. The vicinity of these points are 
the 'valleys', which serve as an additional degree of freedom for excitations with momentum 
close to the Dirac point (low energy excitations): two are associated with each band and 
two are associated with each valley. At 7 = 2, the Dirac points merge at the M— point 
(Fig.lc), and for 7 > 2 a gap forms [11-14]. An effective Hamiltonian is obtained by 
expanding </?fc around the extrema of the bands. Its form depends strongly on the strength 
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FIG. 1: Non deformed (a) and deformed (b) honeycomb lattices, that have two sites in a unit cell 
(white), (c) The first Brillouin zone with the high symmetry points. 
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of the deformation. For 1 < 7 < 2, the effective Hamiltonian is an anisotropic Dirac's 
Hamiltonian [11] 

'Hi = v^p^a^ + VyPyay, (5) 



where Vy = \/3^ta/2, Vx = ta\Jl — 7^/4, and p is the momentum measured from the 
extrema of the bands. As 7 approaches 2, vanishes and one must include high order 
terms in p^. Therefore, at the critical deformation (7 = 2), the Hamiltonian has no linear 
term in px [15] 

H2 = -Vstapyay + \ta?pl.ax. (6) 

Deriving the effective Hamiltonian for stronger deformations that are characterized by 7 > 2, 
we find 

A + t^pl _ _ (Hl^.^p^, (7) 



'H. 



where A = ^(7 — 2). Note that the quadratic term in py may not be neglected compared to 
the linear term, since the dispersion obtained from (7) must coincide with the expansion of 
(4). We also note that when the Dirac points merge, the valley degree of freedom vanishes 
and the number of degrees of freedom is reduced to two. 

In all three cases, the effective Hamiltonian has the general form g{p)(Jx + h{p)ay, where 
g, h are functions of p. Defining F{p) = g{p) — ih{p), the eigenstates are 

x(^Hp) = V^(i ±F/\P\), where \F//3\ = 1, (8) 

and '±' indicate the sign of the propagation constant. In real space, ip^\x,y) = 

Scattering: In order to study the scattering problem, we consider a honeycomb lattice 
with additional refractive index step with a corresponding optical potential. Vat. Eq.(l) 
then transforms: H ^ H — V^. We study the scattering of a wave packet from the second 
band with (3 = —Pq, that is initially located at the region of higher index, and is travel- 
ing towards the interface. The height of the step, Vq, is greater than /5o, mimicking the 
scenario considered by Klein. We consider cases: (i) the step is along the direction of the 
/anjfer hopping parameter (?/— direction), and (ii), the step is in the perpendicular direction 
(x— direction). In both scenarios, we calculate the transmission probability below and above 
the critical deformation, and find it to be qualitatively different. We emphasize that below 
the critical deformation, the system is described by Tii, therefore there is no qualitative 
difference between different directions, i.e., the direction of the step is insignificant. 

Step in y: The additional optical potential is VN{x,y) = Vq ■ Q{y), where 9 is the 
Heaviside function. The transmission (reflection) probability, T (i?), is given by the ratio of 
the transmitted (reflected) current and the incident current: 

R=\Ar\\ T= ,f\A A\\ (9) 
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FIG. 2: (color online) The transmission probability as a function of the angle of incidence, for 
scattering off a step in y (a), and a step in x (b). The various lines correspond to deformations 
above and below the critical value, 7 = 2. 



where At {At) is the transmission (reflection) amphtude obtained from continuity of il) at 
the boundary: 

where qy = [{Vq — /3o)^ + v'^pIY^'^- The negative momentum in X- is due to the fact that 
states from the second band have momentum opposite to their velocity. 

At normal incidence, we find a striking result: the transmission probability is exactly 1 
for all deformations, ranging from a non-deformed lattice and up to the critical deformation. 
Since the waves are transferred to a region with lower index of refraction, and a single 
boundary exist, the result is non-resonant unit transmission in an optical system in analogy 
with Klein tunneling. Note that our result is not sensitive to the wavelength or the height 
of the step, and since there is a single boundary, the unit transmission does not result from 
interference. As such, the unit transmission obtained here is completely non-resonant. This 
is an exceptional case in optics where other cases of unit transmission are resonance effects 
characterized by fine tuning of wavelength and/or potential parameters. 

Another example of non-resonant unit transmission was obtained in graphene (non- 
deformed honeycomb) in Ref. [2] where the unit transmission was linked to the vr geometrical 
(Berry) phase accumulated by circumventing the Dirac point. Our findings do not support 
this linkage, since at the critical deformation the Dirac points merge and the geometrical 
phase vanishes. Nevertheless, the transmission probability at normal incidence is 1, implying 
that it is not directly related to the geometrical phase. Moreover, we find that, for increasing 
7 above 1 the dependence of T on the angle of incidence, 6, decreases significantly, and at the 
critical deformation, T is angle independent up to very large angles (Fig. 2a). Even though 
the transverse wave-vector is conserved, the transverse current changes sign, implying that 
the wave experiences negative refraction. Due to extremely weak angular dependence, the 
reflected wave is negligible, in contrast to other systems where negative refraction is ac- 
companied by significant reflection [17]. At deformations above the critical one, the T{6) is 
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FIG. 3: The input (a) and output (b) intensities for a step in the y— direction, for deformation 
close to the critical one. The position of the step is indicated by the vertical solid line, 
(c) the mean propagation constant along the propagation (dots). The solid line indicates the 
intersection of the bands. 

always smaller than 1, but the angular dependence is still extremely weak (Fig.2a). 

Step in x: The optical potential is VN{x,y) = Vq ■ 9(x). Well below the critical 
deformation, the leading term in the Hamiltonian is linear in p^, and the behavior is identical 
to the previous case, i.e., at normal incidence non-resonant unity transmission is obtained. 
However, as the Dirac points get closer, the Dirac cones are distorted and the valleys are 
no longer separable. This situation occurs when the linear term in px is comparable to the 
quadratic term. As the deformation approaches the critical one, the quadratic term becomes 
the leading term and the nature of the scattering process changes dramatically: since the 
Hamiltonian is quadratic in p^, there are four possible solutions for a specified propagation 
constant (3. Two solutions have a real wave-vector corresponding to propagating waves, 
and the other two have an imaginary wave- vector, corresponding to exponentially decaying 
waves. The presence of the exponential waves changes the transmission probability, and at 
normal incidence the transmission vanishes completely, for all deformations above the critical 
one (Fig.(2b)). This is reflected in a rapid change in T{6 = 0) from unit transmission to 
total reflection over a very small range in the parameters space, e.g., at p^a ~ 0.6, the linear 
and quadratic terms are equal at 7 ~ 1.97, and for higher values of 7 the quadratic term 
dominates. Therefore, we obtain unit transmission for 7 < 1.95, and total reflection for 
7 > 2. The angular dependence is characterized by transmission peaks at small angles, so 
that such potential step filters out the Px = mode. 

Lattices with a deformation close to the critical one, have an effective Hamiltonian that re- 
sembles the Hamiltonian of bilayer graphene, and the scattering properties in the x— direction 
resembles that of bilayer graphene as well [3] . Thus, the critically deformed honeycomb is a 
hybrid of monolayer and bilayer graphene. 
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FIG. 4: The input (a) and output (b) intensities for a step in the deformation direction. The 
position of the step is indicated by the vertical sohd line. 



In order to supplement the analytic treatment which includes various assumptions (e.g., 
tight binding, sharp potential step), we re-examine the scattering problem numerically, using 
a continuous paraxial wave equation with honeycomb photonic lattice [9, 12] (not relying on 
tight binding at all). We solve the eigenvalue problem for deformed honeycomb refractive 
index (Fig. lb) and find the Bloch waves of the system, where the deformation is close to 
the critical one. We construct the initial wave packet from Bloch waves of the second band 
and propagate it using Eq.(l) with the additional smooth step-like optical potential. For 
an index step in the direction, we find that the entire wave packet is transmitted to the 
region of lower index, manifesting non-resonant unity transmission in a 2D system (Fig.3). 
We calculate the mean propagation constant, {H), during the propagation, and verify that 
the wave packet transforms from the second band to the first (Fig. 3c). Moreover, we repeat 
the simulation with potentials of finite width and various shapes, and indeed, the unit 
transmission is independent of the shape and width of the potential. 

As for the x— direction, in order to demonstrate the total reflection, one must use a 
very broad beam in the y— direction, since even waves with very small momentum in y— 
experience significant transmission. We preform such simulations and find that a wave- 
packet that is completely extended in y is totally reflected (Fig. 4). 

In conclusion, we have shown that the scattering of a wave packet in deformed honey- 
comb lattice is extremely unique: below the critical deformation, waves that are normally 
incident upon a potential step in any direction exhibit non-resonant unit transmission, and 
the angular dependence of the transmission becomes very flat as the deformation approach 
the critical one. Moreover, we have demonstrated perfect Klein tunneling with zero geomet- 
rical phase rather that vr. Beyond the critical deformation, the system exhibits non-resonant 
total reflection for waves incident normally upon a step in the deformation direction, which 
indicates that close to the critical deformation the system has some of the unique charac- 
teristics of both monolayer and bilayer graphene. Unlike other systems (graphene or cold 
atoms in optical lattices) photonic lattices in which one can observe the field itself offer us 
a great opportunity to finally have a direct observation of Klein tunneling. 
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